Adiabatic theorems for generators of contracting 

evolutions 



J.E. Avron, M. Fraas, 



Department of Physics, Technion, 32000 Haifa, Israel 



G.M. Graf, P. Grech 



Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland 



June 5, 2012 



Abstract 



We develop an adiabatic theory for generators of contracting evolution on Ba- 
nach spaces. This provides a uniform framework for a host of adiabatic theorems 
ranging from unitary quantum evolutions through quantum evolutions of open sys- 
tems generated by Lindbladians all the way to classically driven stochastic systems. 
In all these cases the adiabatic evolution approximates, to lowest order, the natural 
notion of parallel transport in the manifold of instantaneous stationary states. The 
dynamics in the manifold of instantaneous stationary states and transversal to it 
have distinct characteristics: The former is irreversible and the latter is transient in 
a sense that we explain. Both the gapped and gapless cases are considered. Some 
applications are discussed. 

1 Introduction 

We develop a framework for the adiabatic theory of systems whose evolution is governed 
by a slowly evolving family of linear operators generating a contraction in a Banach space 
[281 [H HH] • More precisely, we study equations of the form 



where L(s) is, for any fixed s, the generator of a contraction semigroup. 

The framework encompasses a wide range of applications from driven stochastic sys- 
tems generated in a Markovian process, through isolated quantum systems undergo- 
ing unitary evolution generated by Hamiltonians, culminating in open quantum systems 
whose evolution is generated by Lindblad operators. 

Adiabatic evolutions have a geometric character. As we shall see, the manifold of 
instantaneous stationary vectors, namely kerL(s), has a distinguished complement, with 
the property that a vector near the former evolves with a velocity in the latter, to leading 




(1) 



order in e. Hence, to lowest order in the adiabatic limit, the vector is parallel transported 
with the manifold. 

Parallel transport may be described more concretely within a particular context, as 
we will show in the next section. For instance, when vectors represent quantum states, 
the instantaneous stationary states are transported like points of a rigid body (see Fig. [I]). 

R(s = 0) 




Figure 1: An example where the set of instantaneous stationary states forms a simplex, 
here a triangle. The extreme points represent the spectral projections Pi{s), i = 0, 1,2. 
Parallel transport rotates the triangle at time s = (triangle whose boundary is the full 
line) to the triangle at time s = 1 (triangle whose boundary is the dashed line) as a rigid 
body. 

We consider both the case where ker L(s) is protected by a gap condition, i.e. is an 
isolated eigenvalue of L(s), and where it is not. 

In the gapped case we give an adiabatic expansion which reveals that the dynamics has 
distinct characters within the evolving subspace of instantaneous stationary states and 
transversal to it. Notably, as we shall see, the motion within ker L(s) is persistent and 
partly even irreversible, whereas the motion transversal to it is transient in the following 
sense: Consider the adiabatic evolution over a finite interval, traversed at a slow rate e; 
assume that the generator is constant near its endpoints and smooth otherwise, and let 
the initial state be stationary. Then the distance of the final state from the manifold 
of stationary states is exceedingly small in e (in fact of infinite order: 0(e N ) for all N), 
whereas the distance covered within the manifold is typically 0(1) and consists in turn of 
two parts: A geometric and potentially reversible part, due to parallel transport, and a 
subleading irreversible correction as large as 0(e) (see Fig. [2J. As we will see by the end 
of the section, this single result entails contrasting physical consequences for isolated and 
open quantum systems. 

In the gapless case we no longer obtain an expansion, however we prove that the 
dynamics of the system is constrained to the manifold of instantaneous stationary states 
and is parallel transported along with the manifold as e — > 0. As an application, it 
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ker L(s) 




ran L(s) 



Figure 2: The straight lines of the upper and lower bundles represent the kernel and the 
range of L(s), as they change with s and attain (left and right) asymptotes when L(s) 
does. The thin curves are the result of parallel transport. The thick curves illustrate the 
motion within the two subspaces, as described in Theorem [6] and Corollary [7} It shows 
the transient nature of the motion in the range: That part is smaller than any power of 
e, when L(s) does not vary. 



generalizes the adiabatic theorem without a gap condition for the Hamiltonian case ( [TU] 
El [35] ) to a class of open quantum systems. 

Although the framework and the theorems are general and independent of the context, 
the geometric interpretation and the implications of the theorems may depend on it. It is 
instructive to illustrate this point for quantum adiabatic theorems. The most familiar ver- 
sion is formulated for the Schrddinger equation, where the state of the system is described 
by a vector in Hilbert space, and x of Eq. (jlj is \ip). An alternate description could have 
been given for the von Neumann equation, where the state of the system is described by 
a density matrix, and x is p, a positive matrix with unit trace. For an isolated system, 
undergoing unitary evolution, the two descriptions are, in principle, equivalent up to the 
loss of an overall phase information in p; nevertheless the elementary formulation and di- 
rect proofs of the "standard" adiabatic theorem tend to refer to the Schrddinger context. 
A formulation in terms of the density matrix is of course a prerequisite towards the for- 
mulation of adiabatic theorems for open quantum systems, where only the von Neumann 
context survives. The unified approach presented here gives such a formulation. When 
applied to the Schrddinger equation, it has a precursor in [8]; when applied to the von 
Neumann equation, in [27] . 

It pays to examine a parallel formulation of the adiabatic theorem for state vectors and 
density matrices in a simple setting. Consider an isolated, finite dimensional quantum 
system whose evolution is generated by a slowly varying self-adjoint and non-degenerate 
Hamiltonian H(s). In the context of pure states, where x = the manifold of instanta- 
neous stationary states are the eigenvectors associated to a distinguished eigenvalue e(s) 
and lie in the kernel of L(s) = i(H(s) — e(s)J. Note that L(s) "knows about" the (instan- 
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taneous) eigenvalue e(s). The adiabatic theorem then says that the evolution within the 
spectral subspace is persistent and depending on history; to lowest order, it is geometric 
and encapsulated in Berry's phase [7] . The evolution transversal to this manifold is tran- 
sient and non-geometric and describes tunneling to eigenvectors of different eigenvalues. 

In the context of density matrices, x = p, the generator of adiabatic evolution L(s) 
is the adjoint action of H(s), namely L(s)p(s) = —i[H(s),p(s)]. This generator, being 
invariant to the replacement H(s) — > H(s) — e(s), has no information on the distinguished 
spectral subspace of H(s). Consequently, the manifold of instantaneous stationary states 
(in the simplest setting we consider) is a simplex whose extreme points are the (instanta- 
neous) spectral projections, see Fig. [TJ In this picture, Berry's phase gets lost; however, 
the associated curvature remains hidden in the motion transversal to the manifold, as 
revealed in some instances of linear response theory, like for the quantum Hall effect (see 
also [5]). 

We finally consider a class of open quantum systems which, though not Hamiltonian, 
preserve the Hamiltonian. The generator of the dynamics, called a dephasing Lindbladian, 
retains the above simplex as its manifold of instantaneous stationary states. If the initial 
data start at a vertex, the motion within the manifold of stationary states simply follows 
the parametrically rotating vertex (see Fig. [T]) — this being the geometric part arising at 
lowest order — ; but to next order the motion is irreversible, non- geometric and directed 
away from the vertex. It is interpreted as tunneling, in the sense of quantum transitions 
between states protected by an energy gap, which may but need not be, a coherent process. 

The anticipated, contrasting consequences are now evident. By the general result 
stated earlier, and under its conditions, for systems undergoing unitary evolution tunnel- 
ing is reversible, since it eventually dwindles to a remainder of infinite order, while for 
systems governed by a dephasing Lindbladian, tunneling is irreversible and comparatively 
large, 0(e). 

The plan of the paper is as follows. In Sect. [2] we describe the general adiabatic 
theorems and the properties of parallel transport. In Sect. [3] we apply these results to 
unitary and Lindbladian quantum systems, as well as to driven stochastic processes. All 
proofs, except for a few short ones, are assembled in Sect. |4j 

2 General results 

In the general scheme mentioned in the Introduction the state space is a Banach space 
and the generators are those of contraction semigroups. We shall present two adiabatic 
theorems which, like their Hamiltonian counterparts, either rely on a spectral gap [HI [19] 
or forgo it [2]. Both depend on the notion of parallel transport. Some preliminaries, like 
the existence of the evolution and the definition of parallel transport, shall be dealt with 
first. 

2.1 Preliminaries 

Propagator. We consider the evolution Q with time-dependent generators L(s), possibly 
unbounded, and state sufficient conditions for the existence of the propagator on a Banach 
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space B. 



Definition 1 Operators L(s), (0 < s < 1) on B are called a C k -family if: L(s) are 
closed operators with a common dense domain D and the function L, taking values in 
the Banach space of bounded operators D — > B, is k-times differentiate in s. Here D is 
endowed with the graph norm of L(s) for any fixed s. 

Lemma 1 Let L(s), (0 < s < 1) be a ^-family and, for each s, the generator of a 
contraction semigroup on B. Then there exist operators U £ (s, s') : B — > B, (0 < s' < s < 
1) with U £ (s, s')D C D, U E (s, s) = 1 and 

e^U £ (s,s')x = L(s)U £ (s J s / )x, (x G D). (2) 



For x G D the unique solution x(s) G D of with initial data x(s') = x is x(s) = 
U £ (s,s')x. Moreover, 

\\U £ (s,s')\\<l (3) 

and 

s-^U E (s, s')x = -U e (s, s')L(s')x, (x G D). (4) 
We will call U £ (s, s')x a solution of ([!]) even for x ^ D. 

Remark 1 By definition ([32], Sect. X.8), a contraction semigroup is strongly continu- 
ous. Its generator is thus closed and densely defined. 

Remark 2 Suppose, in alternative to the hypothesis of the lemma, that the generator 
L(s) is bounded and strongly continuous, and that e = 1. Then the propagator exists and 
is bounded (but not necessarily by 1), uniformly in < s', s < 1 [22]. 

Parallel transport. In adiabatic evolutions the manifold of instantaneous stationary states 
associated to ker L(s) plays a distinguished role. This motivates our interest in families 
of projections. We consider P(s) : B — > B, (0 < s < 1) to be any C 1 -family of projections 
in norm sense. Let P{s) = dP(s)/ds. Then the parallel transport T(s, s') : B — > B is 
defined by 

^-T(s lS ') = [P( S ),P(s)}T(s, S '), (5) 
T( S / ,s / ) = 1. 

It satisfies P(s)T(s, s') = T(s, s')P(s') and hence respects the ranges of P(s). Parallel 
transport is thus a perfect adiabatic evolution: no transitions from the bundle of projec- 
tions P{s) to that of the complementary projections Q(s) = 1 — P(s), nor vice versa. 

A characterization of parallel transport, given in terms of sections x(s) G ranP(s), 
states that the projected velocity vanishes: 

x{s) = T(s,0)x(0) P{s)x{s) = 0, (6) 
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and likewise for Q in place of P. Indeed, for such sections x = Px + Px and Eq. (J5J) 
reduces to 

^-T(s,s')x(s') = P(s)T(s,s')x(s') (7) 



by PPP = 0; hence the contention Eq. (|6j). 

The parallel transport determined by the dual projections P(s)* : B* — > B* is 

T*(s,s') = (T(s',s)y, (8) 

as can be seen from Eq. Observe that unless B and B* coincide, the notion of 

orthogonal projection does not make sense a priori. In the applications both projections, 
orthogonal and otherwise, play a role. 

It is often the case that open systems evolve towards a unique equilibrium state or a 
steady state. This situation is associated with rank 1 projections with special properties 
(see Lemma [4] and Example [3] below) and motivates the interest in this class. 

Lemma 2 Let P{s) be a C 1 -family of rank 1 projections. IfkerP(s) is independent of 
s, then P(s) vanishes on ker P(s) and P{s) = T(s, s')P(s') . 

Note that, without making additional assumptions (e.g. that P(s) is an orthogonal 
projection), parallel transport is not guaranteed to be a contraction. By Remark |2j one 
can only conclude that 

sup ||T(s,s')ll < oo. (9) 

0<s',s<l 



2.2 States 

States of a physical system often enjoy more properties than mere vectors in a Banach 
space. The additional structure we will introduce allows for further geometric properties of 
parallel transport. To mark the difference with the previous and the following subsection, 
we shall denote states by p, rather than by x. The fundamental objects, however, are the 
observables, denoted by a, and their algebra A. 

In the following let B = A* be the dual of a C*-algebra with identity A. We consider 
a second C*-algebra A and bounded linear maps $ : A — > A enjoying 

(i) $ is positive ($ > 0): a > =>■ $a > 0; 

(ii) $ is normalized: $(1) = 1. 

The maps satisfy ||$|| =1 ([H], Cor. 3.2.6) and form a norm closed convex set. For A = C 
one is considering linear functionals, denoted p G ^4*, and (i, ii) define states, p G -4^. 

(The subscripts indicate that the functionals are positive and normalized.) For A = A, 
the dual maps $* : A* — > A* satisfy the corresponding properties (i) p > =3- $*p > 
and (ii) ($*p)(l) = p(l). We call them state preserving maps. By duality, 

11*1 = 1- (10) 
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The maps $ and $* then refer to the Heisenberg and the Schrodinger picture, (with $ 
acting on observables and $* acting on states). We will consider state preserving maps 
which are projections V : A* — > A*. (For economy of notation we omit the star and write 
for the predual, if need arises.) Associated to them are the states in their ranges, 
S := A* + 1 R ran V. Such projections naturally arise through the mean ergodic theorem 
Q16J, Thm. 18.6.1) as projections on stationary states of state preserving semigroups 



POO 

P = lim 7 / e _7 '$,*dt, 



provided the limit exists in norm. 

Remark 3 In case A does not have an identity, we obtain A by adjoining one (\TTj, Def. 
2.1.6). We consider maps defined on A satisfying (i, ii), provided they are compatible with 
the adjunction. More precisely, we consider linear functional p G A* ( and in particular, 
states), provided they arise by canonical extension from p G A (111], p. 52). Of a state 
preserving map it is then required to be so also w.r.t. the amended sense of states. 

Example 1 The compact operators A = Com(7i) on a Hilbert space % form a C* -algebra 
(an identity may be adjoined). Its dual A* = Ji(H) are the trace class operators . Any 
p G j7i("H) is a state if p > and tr p = 1. An example of a state preserving projection V 
is 

v p = J2^p^ en) 

i 

where the Pi are an orthogonal partition of unity on H. As required by the definition of 
state preserving maps, V is the dual of a positive normalized map V* on Com("H). In 



fact, V* also acts by (11) 



The following proposition is concerned with families of projections V(s) and, more 
precisely, with the corresponding parallel transport T(s, s'), determined by Eq. ([5]), and 
states S(s): The action of the former on the latter is that of a rigid motion. In the context 
of Example [T] the proposition is illustrated in Fig. [TJ 

Proposition 3 (Rigid transport) Let the C l -family of projections V(s) : A* — > A*, 
(0 < s < 1) be state preserving. Then T(s, s')V(s') is also state preserving. In particular, 
T(s, s') maps 

• S(s') to S(s) isometrically; 

• (isolated) extreme points of S(s') to corresponding ones of S(s). 

Moreover, if p(s) G ranP(s), depending continuously on s, is an isolated extreme point 
ofS(s), then p(s) = T(s, s')p(s') . 
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Example 2 (continuing Example^. With respect to the Hilbert- Schmidt inner product 
induced by the inclusion J\{%) C J?f(H), the projection V is orthogonal and the transport 
T unitary. These two properties are seen from the following consideration: Since J\CH) C 
Com('H), the action of V* : Com('H) — > Com('H) can be compared with that of V : V* 
preserves J x (n) and V* \ J X (U) = V. We thus have %(s,s') \ Ji(H) = T(s,s') by @, 
besides of (T(s, s')p)(%(s, s')a) = p(a) by 

The example is illustrated in Fig. [T] The motion is rigid in the metrics of both J\(l-L) 
and Jii^C). Explicitly, if p(s) = ^2 \jPj(s) then p(s') = ^jPj( s ') f° r the same Xj, while 
the projections retain their distances in both norms. 

We conclude with a consideration about rank 1 projections, which is linked to Lemma|2} 
In the present setting its hypothesis is satisfied: 

Lemma 4 Consider state preserving projections V of rank 1. Then TdnV* = span{l} 
and kerV is independent of V . In particular, ifV(s) is a C 1 -family of such projections, 
then V(s) = T(s, s')V(s') and p(s) = T(s, s')p(s') , where p(s) is the unique state in 
ranP(s). 

Example 3 Let A = Com('H) and let po £ <Ji(ft) be a state. Then the rank 1 projection 
V : p i-)- (trp)po is state preserving with kerV = {p \ trp = 0}, and V* : a i-> tr(p a)l. // 
Po — po{s) is a ^-family, then V(s)p = (trp)po(s) and the statements of the lemma are 



evident. Note however that, in contrast to the projection (11), the actions of V and V* 
are different. Hence V is not orthogonal in 3?fti\ 



2.3 An adiabatic theorem in presence of a gap 

We assume that is an isolated point of the spectrum of L, which is what we mean by a 
gap. Then, for small e, the differential equation forces a fast time scale of order O^ 1 ) 
on vectors transverse to the null space kerL(s). That scale reflects itself in a fast motion, 
consisting of oscillations and decay. By contrast on vectors in the null space the dynamics 
is slow by x = 0. Nevertheless these vectors leak out of that subspace, because it is itself 
changing with s. The leakage however remains of order 0(e), as shown by Theorem [9] 
below. A complementary result, Theorem[6j constructs a "slow manifold", where solutions 
x(s) remain suitably close to kerL(s) and the time scale is 0(1). Before presenting the 
two results, which are illustrated in Fig. [3j we need to specify the transversal subspace 
complementing kerL(s). 

The general assumptions on L = L(s), (0 < s < 1) are 

Hypothesis 1 L is the generator of a contraction semigroup on a Banach space B. 

As a consequence one has 

Proposition 5 The null space and the range of L, the generator of a contraction semi- 
group, are transversal in the sense that 

kerLHranL = {0}. (12) 

The issue whether the two spaces complement each other is covered by the next hypothesis. 
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Hypothesis 2 The range of L is closed and complementary to the (closed) null space of 
L: 

B = ker L© ran L, (13) 
and the corresponding projections are denoted 1 = P + Q. 

Remark 4 We recall that B\ © B 2 is the notation for the sum B\ + B 2 of subspaces 
Bi C B (i — 1,2) in the case that any vector x in the sum admits a unique decomposition 
x = x\ + x 2 with Xi G Bi. Any two among the statements "Bi (i = 1,2) are closed", 
"B\ + B 2 is closed" , and "Pi : x M- a?j (i — 1,2) are bounded" imply the third. 

Hypothesis 3 L(s) is a C k -family for which remains a uniformly isolated eigenvalue. 

Remark 5 We will see by Hypothesis that zero is either in the resolvent set or an 
isolated point of the spectrum cr(L). In the latter case, by Hypothesis^ the gap is then 
assumed to be uniform. The restriction L \ ranL has a bounded inverse, denoted by L~ L , 
and P(s) and Lis) -1 are C k in norm. 



Remark 6 We will give sufficient conditions for Hypothesis^ in Subsec. \2.5\ For short, 
it is the regular case, given Hypothesis [IJ 

For £ = 0, Eq. Q requires x(s) G ker L(s). For small e the differential equation admits 
solutions which remain close to ker L(s). The construction of the "slow manifold" reduces 
to a differential equation for the slow variables only, with the fast ones providing the 
inhomogeneity. The latter, rather than being governed by a further, coupled differential 
equation, are enslaved to the solution at lower orders. More precisely, the solutions are 
described as follows. 

Theorem 6 (Slow manifold expansion) Let L(s) be a C N+2 -family of operators sat- 
isfying Hypotheses [Ip| Then 



1. The differential equation ex = L(s)x admits solutions of the form 

N 

x(s) = e n {a n {s) + b n {s)) + e N+1 r N (e, s) (14) 

n=0 

with 



• a n (s) G ker L(s), b n (s) G ranL(s). 

• initial data x(0) specified by arbitrary a n (0) G kerL(O), r^(e,0) G B; how- 
ever, the b n (0) are determined below by the a n (0) and together define the "slow 
manifold". 

2. The coefficients are determined recursively through (n = 0, . . . , N) 

bo{s) = 0, 

a n (s) = T(s,0)a n (0)+ [ T(s, s')P(s')b n (s')ds' , (15) 
b n+l (s) = L(sr 1 P(s)a n (s) + L( S )- 1 Q(s)6 n ( S ). (16) 
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3. The remainder is 

r N (e,s) = U E (s,0)r N (e,0) + b N+1 (s) - U E (s,0)b N+1 (0) - / U e (s, s')b N+1 (s')ds 



(17) 

where U e (s,s') is the propagator described in Lemma^ It is uniformly bounded in 
e, ifr N (e,0) is: 



N 



sup \\r N (e,s) || < C N J2\K(fy\\ + IKM)||, 

S n=0 

where Cn depends on the family. 
Explicitly: for a^O) = we have 

ao(s)=T(s,0)a (0), (18) 
b 1 (s)=L(s)- 1 P(s)a (s), (19) 

ai (s)= f' T(s,J)P(s')L(J)- 1 P(s n )a (J)ds'. (20) 
Jo 

Corollary 7 If L(s) is constant on an interval I C [0, 1], then 

b n (s) = 0, (se I). 



Proof. This follows recursively from (16) by P(s) =0. □ 



Corollary 8 If P(s) are rank 1 projections andker P(s) is independent of s, then a n (s) 

tmk(o). 



This is the case in Example |3j See Subsec. 3J2 for an application. 



Proof. In Eq. (15) we have P(s')b n (s') = in view of Lemma [2] and of b n (s') G ranL = 
ker P. □ 



In Theorem [6] the initial data x(0) = P(0)x(0) + Q(0)x(0) is such that the first (slow) 
part is arbitrary, and it prescribes the second (fast) part, up to a remainder. The general 
case that both parts of the initial condition are arbitrary is addressed by a result on the 
decoupling of the slow variables from the fast variables: 

Theorem 9 (Decoupling) Let L(s) be a C 2 -family satisfying the assumptions of The- 
orem^ Then for any solution x(s) of Eq. |7]) 

\\P(s)x(s) -r(s,0)P(0)ac(0)|| < Ce||a;(0)||, (0 < s < 1), 

where C depends on the family. 

Remark 7 No statement about the fast part, Q(s)x(s), is made. The theorem may in 
particular be applied to the difference x(0) = Q(0)x(0) of initial conditions sharing the 
same slow part; in this case, \\P(s)x(s)\\ < Ce\\x(0)\\. 
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Figure 3: The figure shows the result of a computation of the unitary adiabatic evolution 
of a qubit, see Subsec. |3.3| for details. The state is represented as a point on the Bloch 
sphere, Eq. (35). The (red) meridian shows the manifold of instantaneous stationary 
states, i.e. kerL. The parametrization corresponds to uniform speed along this path. 
The "slow manifold" is represented by the (green) curve essentially parallel to the (red) 
meridian. An orbit is shown by the (blue) cycloid. Note that the initial conditions do 
not lie on the slow manifold (61 (0) 7^ when P(0) 7^ 0). This is the reason for the large 
oscillations. 



The proof of Theorem [9] will depend on the following result. We consider linear forms 
<f E B*, the dual of B. The duality bracket is (<p,x). 

Proposition 10 (Adiabatic invariants) Let L(s) be a C l -family as above. Suppose 
the family ip(s) G B* , (0 < s < 1) satisfies 

(p(s) e kei L*(s), ip(s) G ranL*(s). (21) 

Then if is an approximate adiabatic invariant in the sense that for any solution x(t) of 

(<p,x)\ s = e [ (L*- 1 < f ,x)ds'. (22) 

Assuming C 2 -regularity, the expression is bounded as 

Kv.^lSI^Cell^llHxCO)!!, (23) 
where C depends on the family L(s). 

2.4 An adiabatic theorem in absence of a gap 

In the absence of a gap a weaker replacement for the previous theorem is provided by the 
following result, which relies on Hypothesis [T] and replaces Hypotheses [2]{3] by 
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Hypothesis 2' 



B = kerL © ranL. 



(24) 



Hypothesis 3' L(s) is a C 1 -family. 

Theorem 11 (Gapless) Let L(s) satisfy Hypotheses^ and^' for all < s < 1 and 

let P(s), for almost all s, be the projection associated to ker L(s) in the decomposition 



(24); moreover let P(s) be defined for all < s < 1 and C l as a bounded operator on B. 



Then the solution of ex = L(s)x with initial data x(0) = P(0)x(0) satisfies 

sup \\x(s) — r(s, 0)ac(0)|| ->• 0, (e -»■ 0) . (25) 
se[o,i] 

Remark 8 T/ie theorem generalizes the Hamiltonian adiabatic theorem in absence of gap 
fIR UJ/ to the non- self- adjoint case. Actually, the C 2 -regularity of P(s) assumed there is 
relaxed here to C 1 thanks to a remark by Elgart, reported in [35]. 

Remark 9 The "almost all" formulation [10, 35] allows for eigenvalue crossings. 



Proposition 10 has the following variant in the gapless case. 



Proposition 12 Let L(s) be a C x -family satisfying Hypothesis 1. Suppose the family 
<f(s) e B* , (0 < s < 1) satisfies 

cp(s) E kerL*(s), <p(s) = L*(s)<f)(s) (26) 

with uniformly bounded <f)(s) and <p(s). Then 

\(v,x)\ s \<3e snp (\\<f>(s')\\ + ||0(s')||) ||*(0)||. (27) 

0<s'<l 

2.5 Complementarity of subspaces 

In this subsection we will give sufficient conditions for the complementarity Hypothesis [2] 
in relation with a spectral gap, and 2' in its absence. As a help to gauge them, note that 
both are false if is an eigenvalue of L with non- vanishing eigennilpotent, but they hold 
true for a self-adjoint operator L on a Hilbert space if is an isolated resp. non-isolated 
point of its spectrum. 



The two subspaces in Eqs. ( 13 24 ) are transversal, 



kerL D ranL = {0}, 

as a consequence of Hypothesis [TJ as we shall see. However they may fail to generate B 
without further hypotheses. Such hypotheses are given in two lemmas corresponding to 
the two cases. There, a prime indicates a hypothesis tailored to the second, gapless case; 
a sufficient condition for an earlier hypothesis is noted by an added roman numeral. 
Counterexamples matching the two cases are also given. Related results are found in 
Thm. 18.8.3). 



12 



Lemma 13 Let B be a Banach space and L a closed operator on B. Assume, besides 
Hypothesis^ that 

(H2i) If is in the spectrum, o~(L), then is a discrete eigenvalue. 



Then B = ker L © ranL, cf. Eq. (ISfy. 

Property (H2i) implies that L is Fredholm, and hence 

(H2ii) L is semi- Fredholm. 

In conjunction with Hypothesis^ Properties (H2i) and (H2ii) are equivalent. 

Remark 10 By definition, a discrete eigenvalue is an isolated point A of the spectrum, 
the Riesz projection of which, 



— l(L-z)- 1 dz, (2* 

>7T1 J r 



2lT\ 

is finite- dimensional. Here A is the only point of the spectrum encircled by V . 
Remark 11 Hypothesis (H2i) is trivially satisfied if dim. B < oo. 

Remark 12 We recall that L is semi-Fredholm iff ran L is closed and kerL or Bj ranL 
are finite- dimensional. If both are, L is called Fredholm. 



Example 4 Assumption (H2i) can not be omitted from Lemma 13 if the splitting (13) is 
to be ensured. In fact in (\2J$ , Thm. 2.2) an example is given of a non-trivial generator 
L of contraction semigroup with trivial null space, yet with o~(L) = {0}. Hence (H2i) fails 
there. By the equivalence with (H2ii), ranL is not closed, spoiling 

Lemma 14 Let B be a Banach space and L a closed operator on B. Assume, besides 
Hypothesis^ that 

(H2'i) B = ker L + ranL. 



Then B = kerL © ranL, cf. Eq. (24). 



Moreover, if ker L + ranL is closed and B reflexive, then (H2'i) follows from Hypoth- 
esis^ 

Recall that, by definition, B is reflexive if B** = B. 

Example 5 Consider the operator L defined by (Lf)(x) = —xf(x) for f G L°°(0, 1) = B. 
Obviously, L has trivial kernel and (e Lt f)(x) = e~ xt f(x), which makes L the generator of 
a contraction semigroup. However, for 1 = g 6 L°°(0, 1) one has 

\\g-Lf\\ LB o >i, (/a»(o,i)). 



Thus ranL is a proper subspace of L°°(0, 1). In relation with Lemma 14 , the example 
shows that when B is not reflexive (H2'i) does not follow from Hypothesis |i| 

Example 6 As a further, similar example consider the operator L : p i-> — i[H, p] defined 
for p G J\(H) = B, where H is a bounded self-adjoint operator on the Hilbert space H. 
Let H have purely continuous spectrum, so that ker L = {0}. On the other hand, tr Lp = 
for any p G JiiH), because ti Hp = tr pH (J33J. Cor. 3.8). Then tip = extends to 
p G ranL, which is thus a proper subspace of JiiH). 
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3 Applications 



Our results apply to a wide range of driven quantum and classical systems. For quantum 
systems we consider evolutions generated either by a Hamiltonian or a Lindbladian. We 
focus on the special class of "dephasing Lindbladians" which are in some sense inter- 
mediate between Hamiltonians and generic Lindbladians. As we shall explain, adiabatic 
evolutions in the Hamiltonian setting have a different character from those in the dephas- 
ing setting. In the Hamiltonian case tunneling is reversible while in the dephasing one it 
is irreversible. 

In classical systems we consider continuous-time Markov processes. We give an adi- 
abatic expansion for a slowly driven Markov process with unique stationary distribution 
and then restrict our attention to reversible processes and to the generation of probability 
currents. 



3.1 Unitary evolutions 

The results of Sect. [2] may be applied to recover known facts about the unitary adiabatic 
evolution driven by smoothly varying family of self-adjoint Hamiltonians H(s) on a Hilbert 
space % [13 [27j. Consider a simple, discrete eigenvalue e(s) of H(s). Its normalized 
eigenfunction if)(s) spans the manifold of instantaneous stationary states, i.e. the kernel 
of 



L(s) = -i(H(s)-e(s)). (29) 

Eq. ([I]) is the Schrodinger equation and let ip £ (s) be its solution with initial data ip e (0) = 
ip(0). Tunneling, T(s), is defined as the leaking out from the manifold of stationary states, 

T £ (s) = l-ms),Ms))\ 2 . (30) 

There is extensive literature (see [15J and references therein) which is concerned with 
estimates of the tunneling amplitude at all orders in s, or beyond. The simplest version 
of these results can be seen to be a consequence of Corollary [7| Namely: 

Theorem 15 Suppose that H(s) = H(s)* is a C°° -family in the sense of Definition^ 
and is in addition constant near the endpoints s = and s = 1. Let ip £ (s) and ip(s) be as 
above. Then T e (l) = 0(e k ), for any k, see Fig. m 



It may be instructive to examine this result from the perspective of the evolution of 
density matrices generated by the adjoint action of H. In contrast with the Schrodinger 



generator of Eq. (29) whose kernel is one dimensional, the adjoint action —i[H(s),p] has 



a large kernel spanned by all the stationary states. Tunneling is then described not by 



the motion in the range but rather by the motion in the kernel. Theorem 15 may then 
be interpreted as the statement that the adiabatic map is a rigid map of the kernel up to 
terms of infinite order, at time one, see Fig. [T} 

Interestingly, when Theorem [6] is applied to open quantum systems, described by 
a dephasing Lindbladian one reaches the opposite conclusion, namely that tunneling is 
irreversible and 0(e). To complete the picture, it may be worthwhile to discuss from the 
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latter perspective why tunneling gets reduced from first to infinite order in e in the special 
case of the adjoint action. These two points will be addressed in detail in Corollary [19] 
and thereafter. 



3.2 Evolutions generated by Lindblad operators 

Lindbladians arise as generators of dynamical semigroups |23j [13] . Different settings are 
available in the literature. We choose one of them: Let A be a C*-algebra with identity 
and let $t, (t > 0) be a norm-continuous semigroup of positive normalized maps on A. 



As noted in Subsec. 2.2, $ t and have norm 1, and hence are contraction semigroups, 
with dual generators, £* and C. 

Following (23], $ t is called a dynamical semigroup if, besides of the above properties, 
it is completely positive. We find it convenient to follow the accepted tradition and call 
the generator in the Schrodinger picture, £ = (£*)*, the Lindbladian. 

More precisely, the generator in the Heisenberg picture is a (weak-* continuous) 
operator on the Banach space B(%) of bounded operators on a Hilbert space Tl. It is thus 
determined by its restriction on A = Com('H), the compact operators on % (cf. examples 
of Sect. [2]). Then the Lindbladian C : Ji(H) — > Ji(H) has a general form [23] 

Cp = -i[H, P} + \Y1 ([*>> r «] + F*> <<]) ( 31 ) 

a 

with H = H* and r* T a bounded operators on Ti. 

Remark 13 £ does not determine T a , H uniquely. The "gauge transformation" T a — > 
T a + (3 a t and H — > H + i ^2 a {^Q^* a ~ Pa^a)/^ leave C invariant. 

In the generic case the Lindbladian has a 1-dimensional kernel, with ker = span{l} 
independently of C, cf. Lemma [4] We consider a smoothly varying family of Lindbla- 
dians. Let p(s) be the corresponding state and p E (s) be the solution of the adiabatic 
evolution equation Eq. Q with initial data p e (0) = p(0). The tunneling Eq. (30) should 
be generalized to T = 1 — F 2 , where the fidelity is 

F £ (s)=tr{(p(s) 1 / 2 p £ (s)p(s) 1/2 ) 1/2 )- (32) 

In the presence of a gap a system relaxes to its equilibrium state exponentially fast. 
A gapped system with a unique ground state will remain close to the instantaneous 
equilibrium state under adiabatic deformation. For the more interesting, gapless case, see 



Subsec. 3.6 From the results in the previous section we have: 



Theorem 16 Let C(s) be C°° -family of Lindbladians having a unique instantaneous sta- 
tionary state p(s). Then, by Lemma ^ p(s) is parallel transported. Suppose that the 
assumptions of Theorem^hold and is L(s) is constant near the endpoints s = 0, 1. Then, 
by Corollary^ p e (l) = p(l) + 0{e k ) for any k and the tunneling out of the ground state 



%s TJ1) = 0(e k ' 
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3.3 Dephasing Lindbladians 

We say that £ is a dephasing Lindbladian (corresponding to a given Hamiltonian H) if 

ker£* D ker([H, •]) (33) 

as subspaces of B((H). 

By the following proposition, the evolution shares the manifold of stationary states 
with the corresponding Hamiltonian evolution. 



Proposition 17 In connection with Eq. (31) we have: 

1. ker£* D ker([if, •]) is equivalent to Y a = f a (H) for some functions f a . 

2. Y a = f a (H) implies ker£ = kei([H, ■]) as subspaces of J\{H). 

3. If the spectrum of H is pure point, then the last implication is an equivalence. This 
applies in particular to the finite-dimensional case. 

A dephasing Lindbladian conserves all observables which are conserved by H, in particular 
H itself and r a . If one interprets the energy of the system in terms of H and V (see 
Remark 13) then one learns that although the system is open, it does not exchange 
energy with a bath. However, the dephasing Lindbladian induces decoherence w.r.t. the 
energy eigenbasis. A (non-rigorous) scenario where that may arise is discussed in 



Example 7 The simplest dephasing Lindbladian is a 2-level system (a qubit). It is a 
4-parameter family: The Hamiltonian is determined by the 3-vector b 

2H = b-a, (bER 3 , 0= ((71,(72, <7 3 )) 

where o~j are the Pauli matrices and 7 > characterizes the dephasing 

Cp = -i[H, p] + [[H, p),H], (7 > 0). (34) 

(Recall that by 4H 2 = (b ■ b)t any function of H is of the form f(H) = aH + /31; the 
dephasing term is written in such a way that 7 is dimensionless.) The canonical map of 
normalized states into the Bloch ball, 

p H> n E M 3 , \n\ < 1 : p= , (35) 

maps the evolution equation p = Cp into the Bloch equation [1^] 

n = b x n + 76 x (b x n), (36) 

where b = b/\b\ . 
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3.4 Adiabatic expansion for dephasing Lindbladians 



For simplicity consider H with simple eigenvalues eo, 
tors ipi '■ 



, e^-i with normalized eigenvec- 



H 



on a finite dimensional Hilbert space "H, dim'H = d. The operators E^ := l^)^) 
form a basis of <B("H), the linear maps on "H, which is orthonormal once that space is 
endowed with the Hilbert-Schmidt inner product. A straightforward computation using 
Prop. 17 shows that are eigenvectors of £ and the eigenvalues CEij 
< and Xij = if and only if i = 
'i 7^ j) with the corresponding projections (cf. (11)) 



A 



A-,,-, Re A ? ; 



and ran £ by i?^ , 



satisfy 

Hence ker £ is spanned by En = Pi 



X-E 

l 3 l 3 




Figure 4: The states of a qubit (2-level system) can be represented as the 3D ball, the 
interior of the Bloch sphere. For a dephasing Lindbladian, the set of stationary states is the 
(blue) axis whose extreme points (red dots) are spectral projections for the Hamiltonian 
H. In the adiabatic setting the (blue) axis moves slowly. 



We now consider a smooth family of Lindbladians £(s) of dephasing form. 
Theorem 18 The equation 

ep{s) = C(s)p(s) 

admits a solution of the form 
p(s) = P (s) + e ( ^ + ^ J " e £ ( P o( s ) " P M f a M ds ' + ( 37 ) 



with 



-2ReA 0i (s)) 



aj (s) = tr(P (s)P,-( S ) 2 Po(s)) " V , x * • 



More generally, the expansion applies to any solution for which it does at s = 0, e.g. for 
the one with initial condition p(0) = Pq(0), if Pq(0) = 0. 
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The expansion (37) is just p(s) = ao(s) + e(bi(s) + ai(s)) + 0(e 2 ), in this order, with 



coefficients given in (18-20). Like in the Hamiltonian case, &i(s) G ran Q(s) describes 
the shift of the slow manifold relative to the manifold of instantaneous stationary states 
which is reversible in the sense of Corollary [7j Unlike there, ai(s) G ran"P(s) now describes 
irreversible tunneling by means of a loss and a gain term involving Pq(s) and Pj(s), (j ^ 0) 
respectively. More quantitatively, tunneling out of Pq(s) is given by Eq. (J32j) (with p(s) 
there replaced by the rank 1 projection Pq(s)) as 

T £ (s) = l-F 2 (s) = l-tr(p(s)P (s)). 

For arbitrary Pq(0) we have the following result: 

Corollary 19 The solution of ep(s) = C(s)p(s) with the initial condition p(0) = Po{0) 
tunnels like 

ps 

-2\ 



T £ (s) = e V / atjWds' + 0(e 2 
;J „ Jo 



with ctj(s) > 0: Tunneling occurs at a non-negative rate, is irreversible and 0(e). 

This result should be contrasted with the small tunneling of infinite order in the unitary 



case, Theorem 15 Alternatively, that case can be analyzed on the basis of C(s)p = 
—i[H(s),p], following j27j. The solution clearly remains a projection from p(0) = Pq(0) 
on, i.e. p(s) = p(s) 2 . Using the expansion (14) for x(s) = p(s) then yields 



a n + b n = '^2(a j a n - j + ajb n _j + bja n -j + bjb 



3=0 



In view of ao{s) = Po(s), see Example [Tj and bo(s) = 0, this reads a n = c n + a n Po + Poa n , 
where c n depends on Oq, . . . , a n _i, b Q . . . , b n . Since Va n = a n by definition of a n , we obtain 



two recursions for Pia n P{, one for i = and one for i ^ 0. Together with Eq. (16), 
all coefficients a n (s),b n (s) are now determined instantaneously in terms of H(s) and its 
derivatives, and in particular without reference to the history H(s'), (s' < s). As a result, 
the tunneling is of infinite order. 

Example 8 (continuing Example^ The adiabatic expansion Eq. (31) takes a rather 



simple form for the Bloch equation ( p 'by . With h replaced by eh and initial condition 
n(0) = —6(0) one finds 

= - s « + m ( 7(s) ' ( i } + + b(s) + b{s) 1 a(t)it 1 + ° (£% (3f 

where 

lit) |6(t)| 2 



a(t) 



l + 7 2 (t) |6(t)|- 

The terms in parentheses, in the order as they appear, have the following interpretation: 
The first term, being proportional to jb(s) describes friction that causes lagging behind the 
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driver b. The second term describes "geometric magnetism" , a term introduced by JJfl . The 
third term is tunneling and describes motion along the axis towards the center, see Fig. [^} 
While the first two terms describe instantaneous response in the plane perpendicular to 
the stationary axis b(s), the last term describes irreversible motion inside the Bloch sphere 
along the axis, Fig. 0. 



Proof. That (34) defines the most general dephasing Lindbladian follows from its spectral 
properties, since the Lindbladian is uniquely determined by kernel and an off-diagonal 
eigenvalue A, Re A < 0; if stands for the ground state, then ImA > for A = Aoi- The 
Bloch equation follows from the commutation relations [n\ • a, • a] = 2i(ni x nz) ■ a. 
To get the expansion (38) write (37) in the form 

p{s) = P (s) + J-(Bb\{P u P } + ilm\[P 1 , P }) 

- e(P (s) - P^s)) [ a x {s')ds' + 0(e 2 ) 
Jo 

and use A = |6|(i — 7) as well as the ( ant i-) commutation relations 

1; . j . x . \h\ 2 

{Pi, Po} = ~^-a, [Pi, P } = --(bxb). a, (Pi) 2 = ^ 



to get the first order correction terms exactly in the same order as they appear in (38). 

□ 

Solutions of the Bloch equations are illustrated in Fig. [3] 

Further applications of driven dephasing Lindbladians are described in (31 IH [2] . 



3.5 Driven Markov processes 

Theorem [6] may be applied to an evolution of the probability distribution of a continuous- 
time Markov process. In particular, we shall describe below an application to (stochastic) 
molecular pumps [31] (see also [291 [T7] ). 

Let X be a random variable on a finite state space S = (1, 2, . . . , d) and denote 

Pi = Prob(X = i). 

The evolution of X is governed by 

d 

Pi = ^2 L ijPj, ( 39 ) 
3=1 

where the transition rate j —> i, (i ^ j), is non- negative and Ljj := — ^ij- The 
transition matrix := exp(Lt) is a left-stochastic matrix (0 < 4>ij < 1, Y2t=i 'Pij = 1)> 
a contraction in the norm \\p\\i = \P3\1 an d converges to a projection, <p(t) — > P + , 
(t — > 00) ([M], Thm. 4.4.8). The range of P + is spanned by stationary probability 
distributions, meaning Ylj=i ^ J ij 7[ j = 0- 
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We assume that the state space S is indecomposable and denote by it the unique 
stationary distribution of L, whence kerL = span{7r} and ranL = {p | YlPi = 0}- m l me 



with Eq. (13), let P be the rank 1 projection associated to that pair of subspaces, which 
are left invariant by L. We identify L -1 with the map (1 — P)L~ 1 (1 — P) defined on all 
of C d , and denote its matrix elements by L^- . 

Now we consider a smooth family of generators L(s) with corresponding stationary 
states 7r(s). 

Theorem 20 Assume that S is indecomposable for L(s) and that-fr^O) = 0. The solution 
of 

d 

i=i 

with initial condition Pi(0) = 7Tj(0) is 



d 



Pt (s) = 7r,(a) + £^^1(5)^(5) +0(e 2 ). (41) 

3=1 



Proof. The expansion (41) is just that of Theorem [6] Note that by ker P(s) = ranL(s) 
(or, more abstractly, by Lemma [1] for A = £°°(S)) the hypothesis of Corollary [8] is satis- 
fied. Thus T(s, s')tt(s') = 7r(s) and ai(s) = 0. □ 

We say that L satisfies a detailed balance if 

Mr, : l-ij-j (42) 

is a symmetric matrix for some tt, in which case that is the stationary distribution. This 
can be interpreted as the statement that the current through any link j — > i 

Jij{p) = Lijpj - LjiPi (43) 

vanishes at equilibrium, Jij(iv) = 0. 

We now strengthen the assumption on S from indecomposable to irreducible, meaning 
that iTj > 0. Then kerM = span{(l, 1, . . . 1)}, ranM = ranL, and the two subspaces 
decompose M as a linear map. At first M~ l is defined on ranM, and it may be extended 
afterwards, arbitrarily but linearly, to all of C d , e.g. by having it vanish on ker M. 

In applications to (stochastic) molecular pumps one is interested in systems that carry 
no current in their equilibrium states, but can be induced to yield net particle transport 
in an adiabatic pump cycle. Note first that M and tt provide natural coordinates for 
those irreducible processes L which satisfy a detailed balance condition. We set the pump 
period (in scaled time) to be unity. 

The net transport across the link j — y i is expressed in terms of the integrated proba- 
bility current 



1 f 1 

Tij := - Jij(p(s)) ds. 
£ Jo 



The following describes the current in the adiabatic limit. 
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Corollary 21 Let s \-t {M(s),tt(s)} be a pump cycle with tt(s) the unique equilibrium 
state for every s. Assume that 7Tj(0) = 0. Then the transport is geometric to leading 
order, given by 

Ta = f Y, {M^s)M7 k \ s ) - M ri {s)M7 k \ s )) dit k (s) + 0(e). (44) 
Jo k=i 

In particular Ty = 0(e) if tt is constant or, in the periodic case L(0) = L(l), if M is. 
Remark 14 Here, geometric means that the transport is independent of the parametriza- 



tion of the pumping cycle. This is evident in Eq. (44) 



Remark 15 The corollary says that effective pump cycles require the variation of both 7r 
and M. As a matter of fact, the long time average ofTij vanishes under the conditions 
stated in the last line of the corollary regardless of the adiabatic limit f31[ \2h\ ESJ/. 

Proof. The contribution to of order e~ x vanishes due to the detailed balance condition. 



To next order Eqs. (43, 41) yield 



Tij = r$3( L «( a ) L 7*( a ) " L 3l (s)L7 k \s))Tt k (s) ds + 0(e). 
Jo k=i 



Eq. (42) may be written as M = (1 — P)L(1 — P)U, where Uij = ^5^-, implying L 1 



4 - P)UM-\1 - P). Thus, with P jt = tt,-, we have 

LijLjifik = Lij ^(1 - P^M^jCk = M l3 M; k l i, k ~Y,M m M^ k . 
i i 

After interchanging i,j and taking the difference, the second term cancels and we are left 



with (44). The additional claim in the periodic case follows by the fundamental theorem 



of calculus. □ 
3.6 Remarks about the gap less case 



Theorem 11 can be applied to evolutions generated by either a Hamiltonian or a Lindbla- 



dian, just like Theorem [6] was in Subsecs. 34^ andj^, respectively. The Hilbert space H 



must, of course, be infinite-dimensional in order for the gap to vanish. 

In the unitary case the result provides a new proof of the adiabatic theorem for Hamil- 
tonians without spectral gap, as in [21 [35] ; in fact, Hypotheses[l]and[2j are trivially satisfied 
in this case. 



In the Lindbladian case it would be desirable to apply Theorem [TT] to the natural 
space B = Ji(H). Unfortunately its Hypothesis |2j is typically not satisfied, as Example [6 
shows. However, Lindbladians of the dephasing kind have extensions to Jiijl), to which 
that hypothesis does apply. The result so obtained is applicable to tunneling, as we shall 
see below. A typical situation leading to a gapless (dephasing) Lindbladian arises from a 
Hamiltonian having both continuous and point (e.g. discrete) spectrum. 
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Theorem 22 Let C(s), < s < 1, be C 1 -family of dephasing Lindbladians acting on 



J2{H). Then Eq. (24) holds true. Let V(s) be the associated projections, for almost all s; 
moreover let V(s) be defined for all < s < 1 and C 1 as a bounded operator on ^(H)- 
Then the solution of ep = £(s)p with initial data p(0) = P(0)p(0) satisfies 

sup \\p(s)-T(s,Q)p(0)\\ j2m ^Q, (e^O). (45) 
se[o,i] 

The assumptions on the dephasing Lindbladian C(s) follow from corresponding ones on 



the underlying Hamiltonian H(s), see Eq. (33): Let Pj(s) be its eigenprojections, for 



almost all s, cf. Remark |9j and let them be defined for all < s < 1 and C 1 uniformly 
in j. Then, as will be proved together with the theorem, T(s,0)Pj(0) = Pj{s). This 
implies the following result about tunneling out of an initial state given by a rank 1 
eigenprojection p(0) = Po(0): 

T £ (s) = 1 - tr(p(s)P (s)) = 1 - tr((7>, 0)P o (0))P o (s)) + o(l) = o(l). 

The next example illustrates that, although the tunneling out of the initial state is of 
order o(l), an adiabatic invariant is conserved up to order e. 

Example 9 Consider the Hamiltonians H(s) = V(s)HV*(s) arising from a C 2 -family 
of unitaries V(s) and from a bounded H. Let p(s) solve the equation ep = —i[H,p]. Then 
the energy is an adiabatic invariant in the sense that 

|tr(fT(l)p(l))-tr(fr(0)p(0))| = 0(e). 

This follows from Eq. (p??|). We may in fact apply that estimate to x(s) = p(s), (p(s) = 



H(s), (ip, x) = tr(Hp) and C = — i[H, ■], since the assumptions (26) hold true by 



C{s){H{s)) = 0, H{s) = -[H(s), V(s)V*(s)) = C*(s)(iV(s)V*(s)). 

4 Proofs and supplementary results 

We begin by recalling the Hille-Yosida theorem ([32J, Thm. X.47a): A densely defined, 
closed operator L on B generates a contraction semigroup iff 

(0,oo) C^L), 

||(L-7)x|| > ^\\x\\, (7 > 0,x G D(L)). V ' 



Conditions (46) reflect the connection between the resolvent and evolution operators. For 



example the only if part of Hille-Yosida theorem follows from the formula 

poo 

-(L-7)- 1 = / e^-^dt, (7>0). 
Jo 

Proof of Lemma [l| The hypotheses are a convenient strengthening of those of [32] . 
Thm. X.70, including the remark thereafter. All our statements but uniqueness and 
Eq. ^ are among its claims, and those two are consequences of its proof. Alternatively, 
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the results may be read off from [20]: Eq. ^ from Thm. 4, Eq. ^ and uniqueness from 
Thm. 1, and Eq. Q from Thm. 2. □ 

Proof of Lemma [2j Rank 1 projections P are of the form Py = a(y)x where x G B 
and a G B* are determined up to reciprocal factors. Any a with ker a = ker P may thus 
be picked, and then x normalized by a(x) = 1. Since ker P(s) is independent of s, so is 
our choice of a in P(s)y = a(y)x(s), while x(s) is C 1 . Thus P(s)y = oc(y)i(s), which 
vanishes for y G kerP(s). The claim just proved states P = PP; together with (7b both 
sides of P(s) = T(s, s')P(s') are seen to satisfy the same differential equation in s. □ 



Consider the ranges of P(s) and of Q(s) = 1 — P(s). As the name suggests, parallel 
transport T(s, 0) is obtained by projecting vectors from either subspace at to the cor- 
responding one at s or, more precisely, by repeating the procedure on the intervals of an 
ever finer partition of [0, s]. In fact, 

P(s)P(0) + Q(s)Q(0) = 1 + [P(0), P(0)]s + o(s), (s 0), 

implying by Eq. ^ 



N-l 

T(s,s')= lim TT(P( S<+1 )P( Si ) + g(si+i)g(si)) 

i=0 

N N 

hm (n^(*i)+n^))' ( 4? ) 



j=0 i=0 



where s' = sq < s\ < ... < sn = s is a partition of [s',s] into intervals of length 
\s i+ i - s» | = iV _1 |s - s'| and Ililo 1 ^ = • • • Aq. 

Proof of Proposition |3j Since the product of state preserving maps is state preserving 
the first claim follows from (47). As a result T(s, s')V(s') maps S(s') — > S(s), with inverse 



T(s' , s)V(s). The first bullet follows from Eq. (10) for the two maps, i.e. 



\\p(s) ~ ~p(s')\\ = \\T(s',s)V(s)(p(s)-p(s))\\ < \\p(s) - ~p(s)\\. 

Any convex decomposition of T(s, s')p(s'), (p(s') G S(s')) entails one of p(s'), which yields 
the second bullet in the variant where the bracketed word is omitted. The continuity of 
T(s,s')p w.r.t. p yields the other variant. To obtain the last statement we note that 
7~(s, s')p(s') is, for fixed s and all s', an isolated extreme point in S(s), just like p(s). 
They agree, since they do for s' — s. □ 



Proof of Lemma |4j By the setting of Subsec. 2^ V has a predual P*, which is also 



of rank 1. Since = 1 by the normalization condition, we have ranP* = span{l}. 

Thus kerP = (ranP*)- 1 is independent of P. We recall that S L C B = A* is the annihila- 
tor of a subspace S C A. The remaining claims follow from Lemma [2] and Proposition |3jn 

Proof of Proposition [5| La = and a = Lb imply (L — 7) (a + 76) = — 7 2 6 and by (46 ) 

7||a + 7&|| < 7 2 ||6|| for 7 > 0. After dividing by 7 we obtain a = in the limit 7 — > 0. □ 



23 



Before proving Theorem 6] we derive a few further consequences of its assumptions 
Hypothesis |2| includes Eq. (12); moreover the pair of subspaces decomposes L 

Ta.nL) 



The 



restriction L \ ranL is closed and has range ran(L \ ranL) = ranL; by (12) it is one- 
to-one. Thus ^ <r(L r ranL). Together with o[L \ kerL) C {0} we conclude that the 
resolvent set contains a punctured neighborhood of 0, which proves the presence of a gap. 



We maintain that the projection P is given by the Riesz projection, see Eq. (28) 
with A = 0. Calling the latter temporarily P we thus claim Pa = a for a e kerL 
and Pb = for b G ranL. The first statement is evident from (28); for the second it 
suffices, by PL C LP, to show that ranP n ranL = {0}. This in turn follows because 
L \ (rani 3 fl ranL) is a bounded operator with empty spectrum; in fact, it is contained 
in cr(L I" ranP) fl o(L \ ranL) = since the first spectrum is contained in {0}, while the 
second is disjoint from it. Finally, we recall the formula for the inverse of L \ ran(l — P) 
(j2U, Eq. (III.6.23)): 



In particular, P(s) and L(s) 1 are C k in norm. 



Proof of Theorem [6j We insert the right-hand side of (14) as an ansatz into and 
equate orders e n , (n = 0, . . . , N), resp. 0(e N+1 ). We find 



Lb = 0, 
a n + b n = Lb n+ \, 
er N + a N + b N = Lr N . 



(n = N-l) 



(48) 
(49) 



In particular, 6o — 0. Note that Qa = implies Qa + Qa = 0, or Qa = —Qa = Pa. 
Similarly, Pb = —Pb. Applying Q and P to (48) yields 



If b n is known, (51) implies 



Pa n + Qb n — Lb n+ i, 
Pa n - Pb n = 0. 

a n = Qa n + Pa n = Pa n + Pb ri 



(50) 
(51) 



the solution of which is ( 15 ) by Eq. and the Duhamel formula. If a n and b n are known, 
b n+ \ follows from (50) and Lemma 13, provided b n is differentiable (see below). All this 



determines bo, do, &i, • • • , ajv-i? ^>n- We then define aN,bx+i by the same Eqs. ( |15| 16) 
which ensures a N + b^ = Lb N+x . Then (49) reads 



with solution 



er N = Lr N - Lb N+1 



r N (e, s) = U £ {s, 0)r N {e, 0) - e' 1 / U E (s, s')L{s')b N+l {s')ds 



o 

s , d 



U E (s,0)r N (e,0) + J (—U £ (s, s'))b N+1 (s')ds' . 
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An integration by parts yields (17) and the bound on the remainder follows from Hy- 
pothesis [l] through Eq. (pjl). Inspection of the recursion relations shows a n ,b n G C N+2 ~ n , 
which provides the required differentiability. □ 

Proof of Theorem [9| For an adiabatic invariant ip Eq. (23) together with (<p(0),x(0)} = 
(<p(s),T(s,0)x(0)) yields 



Ms),x(s)-T(s,0)x(0))\<Ce\\ l p(s)\\\\x(0)\\. 
The claim follows from ||P(s)x|| = sup{| (<p, x) \ \ if G ker L*(s), \\ip\\ = 1}. 



□ 



Proof of Proposition 10 The assumption (21 ) may be phrased differently. The projec- 
tions P* and Q* are associated to kerL* ©ranL*, because the subspaces (13) decompose 
L. Thus Q*ip = 0, P*ip = just means that ip(s) G ranP*(s) is parallel transported: 
v(s)=T*(s,s'Ms'). 
Eq. (22) follows from 

d 



ds 



(tp,x) = (ip,x) + (tp,x) = (L*L* (p,x) +e 1 ((p,Lx) 
= (L*- 1 cp,Lx) + e- 1 {L*cp,x) = e(L*- 1 (p,x) + 0. 



Integration by parts in (22) gives 



( ¥ »,x>|S = e((0 |i r>|S- / (j>(s'),x(s f ))ds' 



(52) 



where <j>(s') = L*(s') 1 Lp(s'). We observe that ||rr(s')|| < ||x(0)|| by Lemma [IJ and 
(p(s) = P*(s)cp(s) by Q. By (§) we see that ||^(s')||, U(s')\\ and \\<j>(s')\\ are bounded by 

□ 



a constant times \\<p(s) 



, proving Eq. (23). 



We now turn to the case without gap. 



Proof of Theorem Let us first dispose of the "almost all" qualifier in the theorem. 
We observe that, by continuity, L(s)P(s) = holds for all < s < 1. In particular, 
L(s)xo(s) = for x (s) = T(s, 0)x(0). The remainder to be estimated is r(s) = x(s) — 
x (s). By Eq. Q it satisfies the differential equation er(s) = L(s)r(s) — ex (s) with 
solution 



r(s) 



U e (s, s')xo(s')ds' . 



By Eqs. (|6| 24) we have xq(s) G ranL(s) for almost all s. This property is reflected in 
the splitting 

x (s) = L{s){L{s) - 7 ) _1 io(s) - l(L(s) - 7 )" 1 i (s) , (7 > 0) , 



motivated by Eq. (53) below. Together they yield 



r Is) 



U £ (s, s')L(s')(L(s') - i)- 1 x (J)d8' + / U e (s, s'h(L(s') - 7 )- 1 i Q ( S / )^ / , 
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where, by an appropriate choice of 7 > 0, the second integral can be made arbitrarily 
small by means of dominated convergence; in fact, uniformly in e due to ||{7 e (s, s')|| < 1. 
It remains to show that, for fixed 7 > 0, the first integral vanishes with e. To illustrate 
the argument, let us temporarily pretend that z(s) := (L(s) — •y)^ 1 x (s) is C 1 . Since 
ed s iU e (s, s') = —U e (s, s')L(s') an integration by parts yields for that integral 

e f d s ,U E {s, s')z(s')ds' = eU £ (s, *')*(*') ISS - £ f U e (s, s')^- f z(s')ds' , 
Jo Jo " s 

and exhibits the desired property for e — > 0. Finally, we get rid of the additional as- 
sumption by amending the argument as follows. We introduce a mollifier j, (j e Cg°(K), 
J j(x)dx = 1) and set js(x) = S~ 1 j(x/S), (5 > 0); we extend xq continuously outside of 
the interval [0, 1]; and split 

z = (L - 7) _1 (i - js * i ) + (L - 7)~ 1 (j, 5 * x ). 

Since x — js * x — > 0, (6 — > 0) and \\L(L — 7) _1 || is bounded, both uniformly in s, the 
first term contributes arbitrary little to the integral, uniformly in e, if 5 is picked small 
enough. The preliminary argument can now be applied to the second term in place of z. □ 



Proof of Proposition 12 Eq. (52) can be obtained from the present assumptions by 



replacing for L* 1 (p in the previous derivation. □ 
Proof of Lemma |13[ Suppose a closed operator L has as an isolated point in its spec- 



trum, with associated Riesz projection P, see Eq. (28) with A = 0. Then P decomposes L 
with <r(L \ ranP) = {0} and a(L \ ran(l-P)) = a(L)\{0} ([21], Thm. III.6. 17), whence 
ranL D ran(l — P) and kerL C ran P. As a result, (H2i) implies that L is Fredholm 
because both its parts are, and in particular that (H2ii) holds, regardless of Hypothesis 

m 



From now on we assume (HI) and (H2ii), with the former implying Eqs. (46, 12). 
By the stability theorem ([21], Thm. IV.5.31) L — z remains semi- Fredholm for z in a 
complex neighborhood of and the index dimker(L — z) — dim(B/ ran(L — z)) is constant; 
moreover the two dimensions are separately constant in a punctured neighborhood U. By 



(46), they both vanish there, and so does the index at z — 0. This has the following 



implications: First, if G cr(L), then it is isolated. Second, the map kerL — > Bj ranL, 



a 1 — y a -\- ranL is one-to-one by (12) and thus onto by the vanishing index. This proves 



ker L + ranL = £>, completing the proof of Eq. (13). 

Finally in order to prove (H2i), we observe that the Riesz projection is given by P, as 
established at the beginning of this section. Thus it is finite-dimensional because ker L is.D 



Lemma 14 is an immediate consequence of the last two statements of the following lemma. 



Lemma 23 Let L be the generator of a contraction semigroup on B. Then 

bew^L lim7(L-7)" 1 6 = 0, (b E B), (53) 

74,0 



ker L fl ranL = {0} (54) 
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If, in addition, B is reflexive, then kerL + ranL = B. 



Proof. We note that j(L — 7) _1 is uniformly bounded in 7 > by (46). By density, 
it thus suffices to prove the direct implication (53) for b = Lx, for which it follows 
from j(L — ■y)~ 1 Lx = j(x + j(L — 7) _1 x). Conversely, set x 7 = (L — 7) _1 &; then 



6 + 7(L- 7 )- 1 6^6. 



Lx^ 

Next, let b be in the intersection (54): we have (L — 7)6 = —76, and hence b = 
— 7(L — 7)~ 1 6, which vanishes as 7 4- 0. 

To prove the last claim it suffices, by the Hahn-Banach theorem, to show that x* E 
(kerL)- 1 - n (ranL) -1 implies x* = 0. Here S 1 - C B* is the annihilator of a subspace S C B. 
We have (ranL) -1 = kerL* and, in the reflexive case, (kerL) -1 = ranL*. The last equality 
is due to S- 1 - 1 - = S ([21], Eq. III. 1.24). The property (0, 00) C p(L) and the uniform bound 



on 7(L — 7) -1 , (7 > 0) are inherited by L*, and so is the consequence (54). We conclude 
that x* = 0. (As a matter of fact, L* is also densely defined ([21] Thm. III. 5. 29) by reflex- 
ivity, and hence is the generator of a contraction semigroup; in particular (e Lt )* = e L *'.) □ 



rems 



This concludes the proofs of Sect. [2] and we pass to those of Sect. |3j Actually, Theo- 
15 and 16 do not require proof, as they are immediate applications of Sect. |2j while 



the results of Subsec. 3.5 are proven there 



Proof of Proposition 17 To begin we write Eq. (31) as 



Cp = -i[H,p] + ^( 21 > r « - {Pi T « T «})-. 



and hence 



C,a =i[H,a] + ^^2 (21>r a - {a, I^rj) . 



It is evident that T a = f a (H) implies ker£* D ker([L r , ■]) and, through [r*,T a ] = 0, also 
ker£ D ker([Lf, ■]). Together with a detail to be supplied momentarily, the three claims 
thus reduce to the following ones: 

1. ker£, D kei([H, •]) implies T a E {f{H)}" = {f(H) \ f(H) E B{U)}" . 

2. T a = f a (H) implies ker£ C ker([L 7 , •]). 

3. If the spectrum of H is pure point and if ker£ D ker([H, ■]), then r a E {f(H)}". 

Here a prime denotes the commutant. We recall ([U], Thm. 2.4.11) that the bicommutant 
is the strong closure of the original *-algebra. Note that {/(L 7 )} is strongly closed. 
The implications 1.-3. are based on the readily verified identity [23J 

£*(a*a) - a*£*(a) - £*(a*)a = / Ja, T a ]*[a, T a ]. (55) 

Q 

1. Let a E ker([L r , •]) = {f(H)} f . Since that subspace of B(7i) is closed under taking 



adjoints and products, each term in the l.h.s. of (55) vanishes by the assumption of the 
present item, implying T a E {f(H)}". 
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2. Under the assumption, £ acts as under the replacement H — > —H, T a — > T* c 



Since tr£(p) = for p G Ji(H), Eq. (55) implies 



- tr p*£(p) - tr p£(p*) = £ tr[p, ^]*[p, 1^]. 

a 

Thus p G ker£ implies [p, T*] = and, by £(p*) = £(p)*, also [p, r a ] = 0. We conclude 
[H,p] = 0. 

3. By the first assumption we can pick finite-rank projections P n , which are sums 
of eigenprojections of H or of subprojections thereof, such that P n — >1. In particular 
[H,P n } = 0. 

If a G JiiT-L) then £*(a) G j7i("H) and, we claim, tr£*(a) = by our second assump- 
tion. Indeed, it implies £(P„) = and hence 

tr £*(a) = limtr(£„ (a)P n ) = limtr(a£(P n )) = . 

n n 

Let now a G Ji(U) R {/(#)}'. Then £(a) = and tr(£,(a*)a) = tr(a*£(a)) = 0. By tak- 



ing the trace of Eq. (55) we conclude [a, T a ] = 0. The conclusion extends to a G {f(H)}' 
since P„a A a. This proves the claim. □ 



Proof of Theorem |18[ Clearly, 



Vp = ^(PipPi + PipPi) 

(< ^ j: 



We note that 



r(s,s')n(s') = f\(s). 

In fact, the l.h.s. satisfies the differential equation ([5]), viz. 

d 



(56) 
(57) 



ds 



p(s) = V{s)p{s) 



p(s') = Pk(s'), 

and so does the r.h.s., since 

PPfc = PPcP + PjPcP = PfcPfc + PfcPfc = Pk 



follows from (19) and (56): 



The claim now follows from (18 20) with ao(0) = Po(0). Indeed, the middle term of (37) 



C-'VPo = C-'Po = C- 1 PjPoPo + PoPoP, 



^ ^joPjPoPo + Kj'PoPoPj - ^joPjPo + Kj~PoPj- 



(5* 
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For the last term of (37) we compute with (58) 



= 5>i(*Wo), 

(with termwise equality) where we have used PjPfc = —PiPk and tr (P^P^Pj) = tr(P Pj P ) . 
Together with (20) the expansion follows. The generalization follows because of the con- 
traction property of the propagator, Eq. (J3|. □ 

Proof of Corollary |19[ The statement evidently applies to the "slow manifold" solution 
(37), which however does not satisfy p(0) = -Pq(O), as required for the present solution. We 



compare the two by means of Remark [7j For their difference p(s) we have ||p(0)|| = 0(e) 
implying \\V(s)p(s)\\ = 0(s 2 ). □ 



Proof of Theorem 22 As noted at the beginning of Subsec. 3.2| a Lindbladian C is the 
generator of a contraction on J\^H) y while is on Com('H); for a dephasing Lindbladian 
the latter also applies to C itself, as noted in the proof of Proposition 17 Then 
interpolation, C generates a contraction also on J%i^C). 



by 



Given that J~2^H) is reflexive, Eq. (24) follows from Lemma 14 by showing that ker L + 
ranL is closed. We will actually show that the two subspaces are orthogonal w.r.t. to 
the inner product (-,•) of J^(H). In fact, for a dephasing Lindbladian the statement 



ker£ = ker([if, •]) from Prop. 17 also holds true as subspaces of J-iffl), as inspection 
of the proof shows. Thus ker£ C ker by (33). Then La = and b = Cb imply 
(a, b) = (£*a, b) = 0. The theorem follows. 

To prove the statements following the theorem let us note that, for almost all s, 



v(s)p = Y j pAs)pp^)i 



(59) 



where the sum runs over the eigenvalues of H(s). In fact, this follows by the RAGE 
theorem ([121, Thm. 5.8), i.e. 



1 

T 



e im pe -im dt _^ 



(T -»■ oo) 



for p G Com('H). It implies p = PjPPj f° r P £ ker £. The converse is obvious. At this 
point the regularity of V(s) is inherited from that of the Pj(s) 
all s. Finally, T(s, 0)Pj(0) = Pj(s) follows like ([57) . 



and Eq. (59) extends to 

□ 
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